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ABSTRACT 

Despite astrophysical importance of binary star systems, detections are limited to those located in small 
ranges of separations, distances, and masses and thus it is necessary to use a variety of observational tech- 
niques for a complete view of stellar multiplicity across a broad range of physical parameters. In this paper, 
we report the detections and measurements of 2 binaries discovered from observations of microlensing events 
MOA-201 l-BLG-090 and OGLE-201 l-BLG-0417. Determinations of the binary masses are possible by simul- 
taneously measuring the Einstein radius and the lens parallax. The measured masses of the binary components 
are 0.43 M Q and 0.39 M Q for MOA-201 l-BLG-090 and 0.57 M Q and 0.17 M for OGLE-201 l-BLG-0417 
and thus both lens components of MOA-201 l-BLG-090 and one component of OGLE-201 l-BLG-0417 are M 
dwarfs, demonstrating the usefulness of microlensing in detecting binaries composed of low-mass components. 
From modeling of the light curves considering full Keplerian motion of the lens, we also measure the orbital 
parameters of the binaries. The blended light of OGLE-201 l-BLG-0417 comes very likely from the lens it- 
self, making it possible to check the microlensing orbital solution by follow-up radial-velocity observation. 
For both events, the caustic-crossing parts of the light curves, which are critical for determining the physical 
lens parameters, were resolved by high-cadence survey observations and thus it is expected that the number of 
microlensing binaries with measured physical parameters will increase in the future. 
Subject headings: gravitational lensing: micro - binaries: general 



1. INTRODUCTION 

Binary star systems are of astrophysical importance for 
various reasons. First, the y compose an important por- 
tion of stars in the Galaxy dAbt & Levvl fl976t lAbl [l98l 
Duqu ennoy & Mayor! Il99ll) and thus theories about stellar 
formation and evolution should account for the binary na- 
ture of stars. Second, binary stars allow us to directly mea- 
sure the masses of their component stars. The determined 
masses in turn allow other stellar parameters, such as radius 
and density, to be indirectly estimated. These physical param- 
eters help us to understand the processes by which binary stars 
form (Goo dwin et aDl2007t feurgasser et al.l l2007l). In partic- 
ular, the separation and mass of a binary system tell us about 



the amount of angular momentum in the system. Because it is 
a conserved quantity, binaries with measured angular momen- 
tum give us important clues about the conditions under which 
the stars were formed. 

Despite the importance, broad ranges of separations, dis- 
tances, and component masses make it hard to detect and mea- 
sure all binaries. Nearby systems with wide separations may 
be directly resolved using high-resolution imaging, while sys- 
tems with small separations can be detected as eclipsing or 
spectroscopic binaries. However, binaries with intermediate 
separations are difficult to be detected by the conventional 
methods. In addition, it is difficult to detect binaries if they 
are located at large distances or either of the binary compo- 
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nents is faint. As a result, samples are restricted to binaries 
in the solar neighborhood and are not complete down to low- 
mass stars. For a complete view of stellar multiplicity across a 
broad range of physical parameters, therefore, it is necessary 
to use a variety of observational techniques. 

Gravitational microlensing can provide a complementary 
method that can detect and measure binaries that are difficult 
to be detected by other methods. Microlensing occurs when 
an astronomical object is closely aligned with a background 
star. The gravity of the intervening object (lens) causes deflec- 
tion of the light from the background star (source), resulting 
in the brightening of the source star. If the lens is a single 
star, the light curve of the source star brightness is charac- 
terized by smooth rising and fall. However, if the lens is a 
binary, the light curve can be dramatically different, particu- 
larly for caustic-crossing events, which exhibit strong spikes 
in the light curve. Among caustic-crossing binary-lens events, 
those with long time scales are of special importance because 
it is often possible to determine the physical parameters of 
lenses (see more details in section 2). The binary separations 
for which caustic crossings are likely to occur are in the range 
of order AU, for which binaries are difficult to be detected by 
other methods. In addition, due to the nature of the lensing 
phenomenon that occurs regardless of the lens brightness, mi- 
crolensing can provide an important channel to study binaries 
composed of low-mass stars. Furthermore, most microlens- 
ing binaries are located at distances of order kpc and thus mi- 
crolensing can expand the current binary sample throughout 
the Galaxy. 

In this paper, we report the detections and measure- 
ments of 2 binaries discovered from observations of long 
time-scale caustic-crossing binary microlensing events MOA- 
2011-BLG-090 and OGLE-201 l-BLG-0417. In §2, we de- 
scribe the basic physics of binary lensing and the method to 
determine the physical parameters of binary lenses. In §3, we 
describe the choice of sample, observations of the events, and 
data reduction. In §4, we describe the procedure of modeling 
the observed light curves. In §5, we present the results from 
the analysis. We discuss about the findings and conclude in 
§6. 

2. LONG TIME-SCALE CAUSTIC-CROSSING EVENTS 

For a general lensing event, where a single star causes the 
brightening of a background source star, the magnification of 
the source star flux depends only on the projected separation 
between the source and the lens as 



u 2 + 2 



uVu 2 + 4 



(1) 



where the separation u is normalized in units of the angular 
Einstein radius of the lens, 6e- For a uniform change of the 
lens-source separation, the light curve of a single-lens event is 
characterized by a smooth and symmetric shape. The normal- 
ized lens-source separation is related to the lensing parameters 
by 
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where f E represents the time scale for the lens to cross the Ein- 
stein radius (Einstein time scale), to is the time of the closest 
lens-source approach, and uq is the lens-source separation at 
that moment. Among these lensing parameters fo, ?e, and uq, 
the only quantity related to the physical parameters of the lens 



is the Einstein time scale. However, it results from the combi- 
nation of the lens mass, distance, and transverse speed of the 
relative lens-source motion and thus the information about the 
lens from the time scale is highly degenerate. 

When gravitational lensing is caused by a binary, the gravi- 
tational field is asymmetric and the resulting light curves can 
be dramatically different from that of a single lensing event 
(Mao & Paczvhski 1991). The most prominent feature of bi- 
nary lensing that differentiates it from single lensing is a caus- 
tic. A set of caustics form a boundary of an envelope of rays 
as a curve of concentrated light. The gradient of magnifica- 
tion around the caustic is very large. As a result, the light 
curve of an event produced by the crossing of a source star 
over the caustic formed by a binary lens is characterized by 
sharp spikes occurring at the time of caustic crossings. 

Caustic-crossing binary-lens events are useful because it 
is often possible to measure an additional lensing parameter 
appearing in the expression of the Einstein radius. This is 
possible because the caustic-crossing part of the light curve 
appears to be differe nt for events associated with source 
stars of different sizes (Domi nik1[l99HlGaudi & Gouldll 19991 
iGaudi & Petterll2002t[Pejcha & Heyrovskyll20()9l) .~B"v mea- 
suring the deviation caused by this finite-source effect, it is 
possible to measure the source radius in units of the Einstein 
radius, p+ (normalized source radius). Then, combined with 
the information about the source angular size, 9*, the Einstein 
radius is determined as 9% = 9+/p+. The Einstein radius is re- 
lated to the mass, M, and distance to the lens, D^, by 



9 E = (kMtt k] ) 



1/2. 



, 1 1 

7T rel = AU 



(3) 



where k = 4G/(c 2 AU), Ds is the distance to the source, and 
7r re i represents the relative lens-source proper motion. Unlike 
the Einstein time scale, the Einstein radius does not depend 
on the transverse speed of the lens-source motion and thus 
the physical parameters are less degenerate compared to the 
Einstein time scale. 

Among caustic-crossing events, those with long time scales 
are of special interest because it is possible to completely re- 
solve the degeneracy of the lens parameters and thus uniquely 
determine the mass and distance to the lens. This is possible 
because an additional lensing parameter of the lens parallax 
can be measured for these events. The lens parallax is defined 
as the ratio of Earth's orbit, i.e. 1 AU, to the physical Einstein 
radius, = D^9e, projected on the observer plane, i.e. 



7T E = ■ 



Trel 



(4) 



With simultaneous measurements of the Einstein radius and 
the lens parallax, the mass and distance to the lens are 
uniquely determined as 



and 



D,= 



M - 



AU 



K7T E 
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respectively (lGouldl l20Q0). The lens parallax is measured by 
analyzing deviations in lensing light curves caused by the de- 
viation of the relative lens-source motion from a rectilinear 
one due to the change of the observer's position induced by 
the orbital motion of the Earth around the Sun (Gould1 H992l 
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TABLE 1 

Telescopes 



event telescope 

MOA-20 1 1 -BLG-090 MOA MOA-II 1.8m, New Zealand 

OGLE Warsaw 1.3 m, Chile 
^FUN CTIO/SMARTS2 1.3 m, Chile 
ftFUN PEST 0.3 m, Australia 
MiNDSTEp Danish 1.54 m, Chile 
RoboNet FTS 2.0 m, Australia 

OGLE-201 l-BLG-0417 OGLE Warsaw 1.3 m, Chile 

ftFUN CTIO/SMARTS2 1.3 m, Chile 
/jFUN Auckland 0.4 m, New Zealand 
fiFUN FCO 0.36 m, New Zealand 
^FUN Kumeu 0.36 m, New Zealand 
^FUN OPD 0.6 m, Brazil 
PLANET Canopus 1.0 m, Australia 
PLANET SAAO 1.0 m, South Africa 
MiNDSTEp Danish 1.54 m, Chile 
RoboNet FTN 2.0 m, Hawaii 
RoboNet LT 2.0 m, Spain 

lRefsdallfl966l) . This deviation becomes important for long 
time-scale events, which endure for a significant fraction of 
the orbital motion of the Earth. Therefore, the probability of 
measuring the lens parallax is high for long time-scale events. 

3. SAMPLE AND OBSERVATIONS 

We searched for long time-scale caustic-crossing binary 
events among lensing events discovered in the 2011 mi- 
crolensing observation season. We selected events to be ana- 
lyzed based on the following criteria. 

1. The overall light curve was well covered with good 
photometry. 

2. At least one of caustic crossings was well resolved for 
the Einstein radius measurement. 

3. The time scale of an event should be long enough for 
the lens parallax measurement. 

From this search, we found 2 events of MOA-20 1 1 -BLG-090 
and OGLE-201 l-BLG-0417. Besides these events, there exist 
several other long time-scale caustic-crossing events, includ- 
ing MOA-20 1 1 -BLG-034, OGLE-20 1 1 -BLG-0307/MOA- 
2011-BLG-241, and MOA-201 1-BLG-358/OGLE-201 1- 
BLG-1132. We did not include MOA-201 1 -BLG-034 and 
MOA-201 1-BLG-358/OGLE-2011-BLG-1 132 in our analy- 
sis list because the coverage and photometry of the events are 
not good enough to determine the physical lens parameters 
by measuring subtle second-order effects in the lensing light 
curve. For OGLE-20 1 1 -BLG-0307/MOA-20 1 1 -BLG-24 1 , 
the signal of the parallax effect was not strong enough to 
securely measure the physical parameters of the lens. 

The events MOA-201 1 -BLG-090 and OGLE-201 l-BLG- 
0417 were observed by the microlensing experiments that are 
being conducted toward Galactic bulge fields by 6 different 
groups including MOA, OGLE, ^FUN, PLANET, RoboNet, 
and MiNDSTEp. Among them, the MOA and OGLE col- 
laborations are conducting survey observations for which the 
primary goal is to detect a maximum number of lensing events 
by monitoring a large area of sky. The /iFUN, PLANET, 
RoboNet, and MiNDSTEp groups are conducting follow-up 
observations of events detected from survey observations. 
The events were observed by using 12 telescopes located in 
3 different continents in the Southern Hemisphere. In Table 
Q] we list the telescopes used for the observation. 



Reduction of the data was conducted by using photome- 
try codes developed by the individual groups. The OGLE 
and MOA data w e re red uced by photome t ry cod es devel- 
oped by lUdalskil (120031) and iBond et alJ (1200 ll) . respec- 
tively, whic h are based on the D ifference Image Analy- 
sis method (lAlard & Luptoni ri998). The /iFUN data were 
processed using a pi peline based on the DoPHOT software 
(ISchechteretan[T993T) . For PLANET a nd MiNDSTEp data 
a pipeline based on the pySIS software ( Albrow et al.l 120091) 
is use d. For RoboNet data, the DanDIA pipeline (Bramich 
2008) is used. 

To standardize error bars of data estimated from different 
observatories, we re-scaled them so that \ 2 P er degree of free- 
dom becomes unity for the data set of each observatory, where 
X 2 is computed based on the best-fit model. For the data sets 
used for modeling, we eliminate data points with very large 
photometric uncertainties and those lying beyond 3cr from the 
best-fit model. 

We present the light curves of events in Figure Q] and 
12 To be noted is that the overall light curves including 
caustic crossings of both events are well covered by survey 
observations. This demonstrates that the observational ca- 
dence of survey experiments is now high enough to char- 
acterize lensing events based on their own data. The parts 
of light curves with 2455880 < HJD < 2455960 were not 
covered because the Galactic bulge field could not be ob- 
served. Although not included in the selection criteria, both 
events showed a common bump to those involved with caus- 
tic crossings: at HJD~2455655 for MOA-20 1 1 -BLG-090 and 
at HJD~2455800 for OGLE-201 l-BLG-0417. These bumps 
were produced during the a pproach of the sour ce trajectory 
close to a cusp of a caustic. [An & Gouldl (1200 ll) pointed out 
that such triple-peak features help to better measure the lens 
parallax. 

4. MODELING LIGHT CURVES 

In modeling the light curve of each event, we search for a 
solution of lensing parameters that best characterizes the ob- 
served light curve. Describing the basic feature of a binary- 
lens light curve requires 6 parameters including the 3 single- 
lensing parameters to, uo, and t^. The 3 additional binary- 
lensing parameters include the mass ratio between the lens 
components, q, the projected separation in units of the Ein- 
stein radius, s±, and the angle between the source trajectory 
and the binary axis, a (trajectory angle). 

In addition to the basic binary lensing parameters, it is re- 
quired to include additional parameters to precisely describe 
detailed features caused by various second-order effects. The 
first such effect is related to the finite size of the source star. 
This finite-source effect becomes important when the source 
is located at a position where the gradient of magnification 
is very high and thus different parts of the source surface ex- 
perience different amounts of magnification. For binary-lens 
events, this happens when the source approaches or crosses 
the caustic around which the gradient of magnification is very 
high. To describe the light curve variation caused by the finite- 
source effect, it is necessary to include an additional parame- 
ter of the normalized source radius, p+. 

For long time-scale events, such as those analyzed in this 
work, it is required to additionally consider the parallax ef- 
fect. Consideration of the parallax effect in modeling requires 
to include 2 additional parameters tte.n, and tte.e, which rep- 
resent the two components of the lens parallax vector tte pro- 
jected on the sky along the north and east equatorial coor- 
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FIG. 1 .— Light curve of MOA-201 l-BLG-090. The notation in the paren- 
theses after each observatory denotes the passband of observation, where N 
denotes that no filter was used. The insets in the upper panel show the en- 
largement of the caustic crossing parts of the light curve. The lower 3 panels 
show the residuals from the 3 tested models. See section 4 for more details 
about the individual models. 

dinates, respectively. The direction of the lens parallax vec- 
tor corresponds to the relative lens-source relative motion in 
the frame of the Earth at a specific time of the event. The 
size of the vector corresponds to the ratio of the Earth's orbit 
to the Einstein radius projected on the observer's plane, i.e. 
7Te = AU/|/e.Ds/(Z)s —Di)], where = Dl#e is the physical 
size of the Einstein radius. 

Another effect to be considered for long time-scale events 
is the orbital motion of the lens (Dominik 1998; lloka et al.l 
1999tlAlbrow et al.ll2000tlPenny et al.ll201 lMShin et al.ll201 ll 
Skowron et al.ll2.01lh . The lens orbital motion affects lensing 
light curves in two different ways. First, it causes to change 
the binary separation and thus the magnification pattern. Sec- 
ond, it also causes the binary axis to rotate with respect to the 
source trajectory. In order to fully account for the lens orbital 
motion, 4 additional parameters are needed. The first two of 
these parameters are ds±/dt and da/dt, which represent the 
change rates of the projected binary separation and the trajec- 
tory angle, respectively. The other two orbital parameters are 
sii and dsn/dt, where s\\ represents the line-of-sight separa- 
tion between the binary components in units of 6% and dsn/dt 
represents its rate of change. For a full description of the or- 
bital lensing paramete rs, see the summary in the Appendix of 
ISkowron et al.l (1201 ll) . The deviation in a lensing light curve 
affected by the orbital effect is smooth and long lasting and 
thus is similar to the deviation induced by the parallax ef- 
fect. This implies that if the orbital effect is not considered, 
the measured lens parallax and the resulting lens parameters 
might be erroneous. Therefore, considering the orbital effect 
is important not only for constraining the orbital motion of 
the lens but also for precisely determining the physical pa- 
rameters of the lens. 

With all these parameters, we test three different models. 
In the first model, we fit the light curve with standard binary 




5820 5822 5824 



5780 5800 5820 5840 



5800 5900 
HJD-2450000 



FIG. 2.— Light curve of OGLE-201 l-BLG-0417. Notations are same as in 
Fig 1. 



lensing parameters considering the finite-source effect (stan- 
dard model). In the second model, we additionally consider 
the parallax effect (parallax model). Finally, we take the or- 
bital effect into consideration as well (orbital model). When 
the source trajectory is a straight line, the two light curves 
resulting from source trajectories with positive (+«o) and neg- 
ative (-Mo) impact parameters are identical due to the symme- 
try of the magnification pattern with respect to the binary axis. 
When either the parallax or the orbital effect is considered, on 
the other hand, the source trajectory deviates from a straight 
line and thus the light curves with +uq and -uq are different 
from each other. We, therefore, consider both the +mq and 
-Mo cases for each of the models considering the parallax and 
orbital effects. 

In modeling, the best-fit solution is obtained by minimizing 
X 2 in the parameter space. We conduct this in 3 stages. In the 
first stage, grid searches are conducted over the space of a sub- 
set of parameters and the remai ning parameters ar e searched 
by using a downhill approach (Dong et al. 2006). We then 
identify local minima in the grid-parameter space by inspect- 
ing the x 2 distribution. In the second stage, we investigate 
the individual local minima by allowing the grid parameters 
to vary and find the exact location of each local minimum. 
In the final stage, we choose the best-fit solution by compar- 
ing x 2 values of the individual local minima. This multiple 
stage procedure is needed for thorough investigation of pos- 
sible degeneracy of solutions. We choose of s±, q, and a as 
the grid parameters because they are related to the light curve 
features in a complex way such that a small change in the 
values of the parameters can lead to dramatic changes in the 
resulting light curve. On the other hand, the other parameters 
are more directly related to the light curve features and thus 
they are searched for by using a downhill approach. For the 
X 2 minimization in the downhill approach, we use the Markov 
Chain Monte Carlo (MCMC) method. Once a solution of the 
parameters is found, we estimate the uncertainties of the in- 
dividual parameters based on the chain of solutions obtained 
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TABLE 2 
Limb-darkening coefficients 



MOA-201 l-BLG-090 
giant clump • 



quantity 


MOA-2011-BLG-090 


OGLE-2011-BLG-0417 


rv 


0.52 


0.71 


r« 


0.45 


0.61 


r, 


0.37 


0.51 


source type 


FV 


Kill 


Teff (K) 


6650 


4660 


Vturb (km s ') 


2 


2 


logg (cm s~ 2 ) 


4.5 


2.5 



from MCMC runs. 



MOA-2011-BLG-090 



Y ^=5655.0 



OGLE-2011-BLG-0417 

^ = 5800.0 
^ ^=5822.0 





FIG . 3 . — The caustic geometries and source trajectories of the best-fit mod- 
els for MOA-2011-BLG-090 (left panel) and OGLE-201 l-BLG-0417 (right 
panel). The small closed figures composed of concave curves represent the 
caustics. The small open circles represent the lens positions. We note that the 
lens positions and the resulting caustic vary in time due to the orbital motion 
of the lens. We mark 2 sets of lens positions and caustics at 2 different times 
of fi and ti. We also note that the source trajectory, the curve with an ar- 
row, is curved due to the combination of the parallax and orbital effects. The 
source trajectory is presented so that the binary axis is set along the horizon- 
tal axis. The dashed circles represent the Einstein rings corresponding to the 
total mass of the binary lenses. 

To compute lensing magnifications affected by the 
finite-source effect , we use the ray-shoo t ing method. 
( Schn eider & Weissl Il986t iKavser et all 119861: IWambsganssI 
119971) . In this method, rays are uniformly shot from the im- 
age plane, bent according to the lens equation, and land on 
the source plane. Then, a finite magnification is computed 
by comparing the number densities of rays on the image and 
source planes. Precise computation of finite magnifications by 
using this numerical technique requires a large number of rays 
and thus demands heavy computation. To minimize compu- 
tation, we limit finite-magnification computation by using the 
ray-shooting method only when the lens is very close to caus- 
tics. In the adja cent region, we use an analytic hexadeca pole 
approximation (Peic ha & Heyrovskyll2009t iGould 20081). In 
the region with large enough distances from caustics, we use 
a point-source magnifications. 

In the finite magnification computation, we consider the 
variation of the magnification caused by the limb-darkening 
of the source star's surface. We model the surface brightness 
profile of a source star as 



Sx ■■ 



Fx 
tt6L 2 



i-r A ( i--cos^ 



(7) 



where Tx is the linear limb-darkening coefficients, Fx is the 
source star flux, and ip is the angle between the normal to the 
source star's surface and the line of sight toward the star. The 




OGLE-201 l-BLG-0417 
giant clump 

source 




FIG. 4. — OGLE color-magnitude diagrams of stars in the fields where the 
lensing events MOA-201 l-BLG-090 (left panel) and OGLE-201 l-BLG-0417 
(right panel) occurred. The red and blue dots represent the centroid of the 
giant clump and the location of the source star for each event, respectively. 
For OGLE-201 l-BLG-0417, we also mark the position of the blended star. 

limb-darkening coefficients are set based the source type that 
is determined on the basis of the color and magnitude of the 
source. In Table [2] we present the used limb-darkening co- 
efficients, the corresponding source types, and the measured 
de-reddened color along with the assumed values of the ef- 
fective temperature, T e s, the surface turbulence velocity, Vturb, 
and the surface gravity, logg. For both events, we assume a 
solar metallicity. 

5. RESULTS 

In Table [3] we present the solutions of parameters for the 
tested models. The best-fit model light curves are drawn on 
the top of the observed light curves in Figures Q] and |2] In 
Figure [3] we present the geometry of the lens systems where 
the source trajectory with respect to the caustic and the loca- 
tions of the lens components are marked. We note that the 
source trajectories are curved due to the combination of the 
parallax and orbital effects. We also note that the positions of 
the lens components and the corresponding caustics change 
in time due to the orbital motion and thus we present caustics 
at 2 different moments that are marked in Figure [3] These 
moments correspond to those of characteristic features on the 
light curve such as the peak involved during a cusp approach 
or a caustic crossing. To better show the differences in the 
fit between different models, we present the residuals of data 
from the best-fits of the individual models. For a close-up 
view of the caustic-crossing parts of the light curves, we also 
present enlargement of the light curve. 

For both events, the parallax and orbital effects are detected 
with significant statistical confidence levels. It is found that 
inclusion of the second-order effects of the parallax and or- 
bital motions improves the fits with A\ 2 = 571 and 2680 for 
MOA-201 l-BLG-090 and OGLE-201 l-BLG-0417, respec- 
tively. To be noted is that the orbital effect is considerable 
for OGLE-201 l-BLG-0417 and thus the difference between 
the values of the lens parallax measured with (tte = 0.40) and 
without (tte = 0.17) considering the orbital effect is substan- 
tial. Since the lens parallax is directly related to the physical 
parameters of the lens, this implies that considering the or- 
bital motion of the binary lens is important for the accurate 
measurement of the lens parallax and thus the physical pa- 
rameters. 

The finite-source effect is also clearly detected and the nor- 
malized source radii are precisely measured for both events . 
To obtain the Einstein radius from the measured normalized 
source radius, p^, additional information about the source 
star is needed. We obtain this information by first locat- 
ing the source star on the color-magnitude diagram of stars 
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TABLE 3 
Lensing Parameters 



parameters 




MOA-20 1 1 -BLG-090 






OGLE-201 l-BLG-0417 






model 






model 






standard 


parallax 


orbital+parallax 


standard 


parallax 


orb ital+parallax 


X 2 /dof 


5207/5164 


4718/5162 


4636/5158 


4415/2627 


2391/2625 


1735/2621 


to (HJD') 


5688.331±0.121 


5691.563±0.187 


5690.409±0.110 


5817.302±0.018 


5815.867±0.030 


5813.306±0.059 


"0 


0.3307±0.0038 


-O.0613±0.0008 


-0.0785±0.0005 


0.1125±0.0001 


-0.0971 ±0.0003 


-0.0992±0.0005 


(days) 


CM 1 0-1-0 7 1 


97Q cfi^O 97 


990 zin^n 9 1 


fiC\ 7d-t-0 OS 




Q9 9A^n 17 


s± 


0.981 ±0.002 


0.536±0.002 


0.606±0.001 


0.601±0.001 


0.574±0.001 


0.577±0.001 


1 


0.611±0.005 


1.108±0.026 


0.892±0.014 


0.402±0.002 


0.287±0.002 


0.292±0.002 


a (rad) 


-0.181 ±0.004 


0.373±0.005 


0.317±0.006 


1.030±0.002 


-0.951±0.002 


-0.850±0.004 


p* (io- 3 ) 


2.89±0.03 


0.54±0.01 


0.78±0.01 


3.17±0.01 


2.38±0.02 


2.29±0.02 






0.205±0.003 


0.159±0.003 




0.125±0.004 


0.375±0.015 






-0.071±0.005 


-0.022±0.004 




-0.111±0.005 


-0.133±0.003 


ds±/dt (yr -1 ) 






-0.031 ±0.007 






1.314±0.023 


rfa/A (yr -1 ) 






1.066±0.005 






1.168±0.076 


£11 






0.137±0.008 






0.467±0.020 


ds^/dt (yr -1 ) 






-0.784±0.008 






-0.192±0.036 



Note. 



HJD'=HJD-2450000. 





TABLE 4 






Physical Lens Parameters 




parameter 


MOA-20 1 1 -BLG-090 OGLE-20 1 1 -BLG-04 1 7 


M, ot (Mq) 


().82±0.02 


0.74±0.03 


M, (Mq) 


0.43±0.01 


0.57±0.02 


M 2 (Mq) 


0.39±0.01 


0.17±0.01 


Be (mas) 


1.06±0.01 


2.44±0.02 


fj, (mas yr -1 ) 


1.76±0.02 


9.66±0.07 


D L (kpc) 


3.26±0.05 


0.89±0.03 


a(AU) 


1.79±0.02 


1.15±0.04 


P(yr) 


2.65±0.04 


1.44±0.06 


e 


0.28±0.01 


0.68±0.02 


i (deg) 


129.43±0.33 


116.95±1.04 


Note. — 


Mtot : total mass of the binary, M| and Mi : 


masses of the bi- 


nary components, Be'- angular Einstein radius, fi: relative lens-source proper 


motion, De'. 


distance to the lens, a: semi-major axis, P: 


orbital period, e: 


eccentricity, i 


'•: inclination of the orbit. 





in the field and then calibrating the source brightness and 
color by using the centroid of the giant clump as a refer- 
ence under the assumption that the source and clump gi- 
ants experience t he same amount of extinction and reddening 
(lYoo etalJ [2004). The measured V/I colors are then trans- 
lated into V/K color by using the V /I—V /K relations of 
IBessell & Brett! (1 1988b and the angular source radius is ob- 
tai ned by using the V/K color and the angular radius given 
by Kervell a et al.l (12004). In Figure |H we present the color- 
magnitudes of field stars constructed based on OGLE data and 
the locations of the source star. We find that the source star 
is a F-type main-sequence star with a de-reddened color of 
(V-I)o = 0.68 for MOA-20 11 -BLG-090 and a K-type giant 
with (V-/) = 0.98 for OGLE-201 l-BLG-0417. Here we as- 
sume that the de-reddened color and absolute magnitude of 
the giant clump centroid are (V—T)o, c = 1 06 and M/ c = -0.23 
dStanek & Garnavichl 19981) . respectively. The mean distances 
to clump stars of -7200 pc for MOA-20 1 1 -BLG-090 and 
-7900 pc for OGLE-2 01 l-BLG-0417 are es timated based on 
the Galactic model of lHan & Gouldl ( 119951) . The measured 
Einstein radii of the individual events are presented in Table 
[4] Also presented are the relative lens-source proper motions 
as determined by fi = Ob/Ie- 

With the measured lens parallax and the Einstein radius, 
the mass and distance to the lens of each event are determined 
by using the relations (0 and ©. The measured masses of 



the binary components are 0.43 M and 0.39 M for MOA- 
2011-BLG-090 and 0.57 M Q and 0.17 M for OGLE-2011- 
BLG-0417. It is to be noted that both lens components of 
MOA-20 11 -BLG-090 and one component of OGLE-201 1- 
BLG-0417 are M dwarfs which are difficult to be detected 
by other conventional methods due to their faintness. It is 
found that the lenses are located at distances De — 3.3 kpc 
and 0.9 kpc for MOA-201 l-BLG-090 and OGLE-201 l-BLG- 
0417, respectively. 

Since full Keplerian motion of the binary lens is considered 
in our modeling, we also determine the orbital parameters of 
the semi-major axis a, period P, eccentricity e, and inclination 
i. We find that the binary lens components of MOA-201 1- 
BLG-090 are orbiting each other with a semi-major axis of 
— 1.8 AU and an orbital period of — 2.7 yrs. For OGLE- 
201 l-BLG-0417, the semi-major axis and the orbital period 
of the binary lens are — 1.2 AU and — 1 .4 yrs, respectively. 
In Figure [5] we present the distributions of the physical and 
orbital parameters constructed based on the MCMC chains. In 
Table|4] we summarize the measured parameters of the binary 
lenses for both events. We note that the uncertainties of the 
parameters are based on the standard deviations of the MCMC 
distributions. 

We note that the blended light of OGLE-201 l-BLG-0417 
comes very likely from the lens itself, implying that the lens 
can be directly observed. Based on the measured mass of 
0.57 M Q , the primary of the binary lens corresponds to a late 
K-type main sequence star with an absolute magnitude and 
a de-reddened color of Mj — 6.0 and (V — 7)o ~ L5, respec- 
tively. Considering the distance of 0.89 kpc and assuming an 
extinction of A/ — 0.5 and the color index of E(V -I) — 0.3, 
the apparent brightness and color of the lens correspond to 
/ ~ 16.2 and (V-/) — 1.8, respectively. These values match 
very well with the location of the blend marked on the right 
panel of Figure [4] implying that the blend is very likely the 
lens. The visibility of the lens signifies this event because it 
is possible to check the microlensing orbital solution by spec- 
troscopic radial-velocity observation. 

6. DISCUSSION AND CONCLUSION 

We reported detections and measurements of 2 bina- 
ries discovered from observations of microlensing events 
MOA-201 l-BLG-090 and OGLE-201 l-BLG-0417. The light 
curves of the events have common characteristics of long du- 
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MOA-2011-BLG-090 



OGLE-201 l-BLG-0417 




0.75 6 0.85 




3.2 3.4 3.6 
D t (kpo) 







0.75 0.8 0.85 0.9 0.95 1 





Studies of M dwarfs are important not only because they 
are the most abundant stars in the Milky Way but also they 
form a link between solar-type stars and brown dwarfs; two 
mass regimes that might exhibit very different multiplicity 
characteristics. Precise knowledge of multiplicity character- 
istics and how they change in this transitional mass region 
provi de constraints on low-mass star and brown dwarf forma- 
tion dGoodwin e t al. 2007: lBurgasser et al.ll2007l) . Despite the 
importance of M-dwarf binaries, only a few measurements 
of the binary frac tion and distribution o f low-mass stars have 
been made, e.g., iDelfosse et alJ ((2004), and the samples are 
restricted to only binaries in the solar neighborhood. As a re- 
sult, there are still large uncertainties about their basic physi- 
cal properties as well as their formation environment. Consid- 
ering the rapid improvement of lensing surveys both in equip- 
ment and strategy, it is expected that the number of microlens- 
ing binaries with measured physical parameters will increase 
rapidly. This will contribute to the complete view of stellar 
multiplicity across a wide range of binary parameters. 




22 0.24 0.26 0.28 0.3 0.32 




12B 129 130 131 
i (deg) 





1.65 1.7 1.75 1.8 1 65 
/a (mas/yr) 




65 0.7 0.75 




112 114 116 118 120 122 
i (deg) 





FIG. 5. — Distributions of the physical and orbital parameters of the binary 
lens systems for MOA-201 l-BLG-090 (left panels) and OGLE-201 l-BLG- 
0417 (right panels). 



rations with caustic-crossing features, which enabled to de- 
termine the physical parameters of the lenses. It was found 
that both lens components of MOA-201 l-BLG-090 and one 
component of OGLE-201 l-BLG-0417 are M dwarfs. There- 
fore, the discovered microlensing binaries demonstrate the 
usefulness of gravitational lensing in detecting and charac- 
terizing binaries composed of low-mass stars. By considering 
full Keplerian binary motion, we also determined the orbital 
parameters of the binaries. For OGLE-201 l-BLG-0417, the 
blended light comes very likely from the lens itself and thus it 
would be possible to check the orbital solution from follow-up 
radial-velocity observation. 
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